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Abstract 

Decay of B s meson into the dileptonic channels of e + e~ and 
are suppressed in the QCD corrected standard model due to the vector 
and axial nature of the coupling to lepton bilinears.We show that 
in SUSY theories with large tan/3, contributions to the amplitude 
arising out of exchange of neutral Higgs Bosons considerably enhances 
the decay of B s into In some region of parameter space, the 

enhancement is more than two orders of magnitude, bringing it well 
within experimental possibilities in the near future. 
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Flavor changing neutral current (FCNC) decays of the neutral B-mesons 
, B° and B° s provide unique testing grounds of the Standard Model (SM) 
improved by QCD-corrections via the operator product expansions Dur- 
ing the last few years, considerable theoretical attention has therefore been 
focussed on decays like B — > K* 7, B — ► X s 7, B — > X^ l + l~ , B s — > 7 7 
and _B S — > 7/+ i~ in view of the planned experiments at B-factories, which are 
likely to measure branching fractions as low as 10 _8 [Q. The calculations of 
the branching ratios of FCNC B-decay processes are sensitive to new physics 
beyond SM.In particular calculations based on Minimal Supersymmetric ex- 
tension of the Standard Model (MSSM) for many of FCNC processes have 
been done Nothing spectacularly different from SM happens to the the- 
oretical predictions but sizable changes in the branching ratios can occur for 
some regions of the parameter space in MSSM 

The MSSM and in general any two Higgs-doublet model has in addition to 
the unknown masses of the new particle content over and above SM, the ratio 
of the vacuum expectation values of the two neutral Higgs fields, tanf3 = 
V2/V2, as parameters From available data , on can only constrain the 
parameter space within certain regions and typically tan/3 can be constrained 
to be tan/3 < 50, i.e. it can be large [0,0. 

In some recent works IJ, it has been pointed out that the processes B — > 
X s r + r~ and B — > X s ji + \i~ provide unique opportunities to distinguish 
SUSY models with large tan/3. Exchange of Neutral Higgs Boson(NHB) gives 
rise to new amplitude which were negligible in SM but which for large 
tanfi can enhance these processes by as much as 200%. The purpose of 
this note is to show that in SUSY-models with large tan/3 NHB exchange 
contributions causes a sizable enhancement of processes like B s — > u + u~ ', as 
much as by two orders of magnitude in some allowed region of parameter 
space. 

Let us start by recalling the result for B s — > u + in QCD-improved 
Standard Model [BJ. The effective Hamiltonian describing this process is : 



n-i aGp \ 
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with A = V t b V t * s ,Pl,r — \ (1 T 75) an d P = P+ + P-, the sum of the 
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momentum of /i + and pT . C 7 (m b ), Cg^(mb) and Cio(ra&) are the Wilson 
Co-efficients, whose values are given by Misiak ||. Since we are considering 
B s , the matrix element of Tl e ff is to be taken between vacuum and \B®) 
state. This can be expressed in terms of the decay constant f Ba @. 

(0 | * 7" 7* H*?> = -fB^B (2) 

(0\s l5 b\B° s ) = -f B m Ba (3) 

and 

(0| s P R b\B° s ) = (4) 

Since p B = p+ + p^L, the Cg term in eqn(|l|) gives zero on contraction with 
the lepton bilinear, C 7 gives zero by (Q) and the remaining C\q term gets a 
factor of 2 mi. 



Thus 



\Keff\B s ) = ^ 7 =x\f Bs C l0 p' B 



nG *\fB.(Cwm) (Hi) (5) 



Decay rate of B s ^ a + ji is 

- / + /-) = a2 ^% w i w:i 2 ca. ( 6 ) 

The presence of the factor mf makes (§) almost vanish for e + e _ and give rise 
to a small Branching ratio of 2.6 x 10~ 9 for fi + fi~. 

Consider now the MSSM. / + l~ production via the neutral Higgs bosons 
give rise to additional amplitudes and these become important for large 
tan/3 0.The neutral Higgs bosons coupling to b-quark and l + proportional 
to mbtanfi and mitanfi respectively for large tanj3 [§]]. As pointed out in 
H mitan/3 can be large ~ mj, or more. In addition, as has been pointed 
out in ||, the b —>■ s transition via chargino-stop loop gives rise to another 
tanf3 factor, so that the b — > s l + l~ amplitude goes like mbmi(tanfl) 3 for 
large tan[3. Thus these additional amplitude effectively do not suffer from 
the mi suppression of expression (|6]).The additional diagrams can be taken 
care of through and additional effective Hamiltonianp] 



H MB = ^£ X[Cqi (g p R h) { J /} + (s - p R h) { J ^ t) ] (7) 
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The co-efficients Cq 1 ,Cq 2 are Wilson co-efficients. There values at M w scale 
can be calculated by evaluating the Feynman diagrams of NHBs contribu- 
tions. The result is || 

C Ql (m w ) = - T bnifl n n tan 2 (3{(sin 2 a + hcos 2 a)[—(fi(x Ht ) - fi{x wt )) 
ttm^psm ■ uw x wt 



i 
i 



^ m w cosfJ ^ 



+ 2X;(s 1 (i,i')r 1 (<,0 +A 1 (i,i')r 2 (i ) i / ))} (8) 

a' 

C Q2 (m w ) = — ™ hm % Q tan 2 p{\-^—{f 1 {x Ht ) - h{x W t)) + 2f 2 {x H t,xwt) 
rnn AO sin l vw xwt 

m w cos/3 

+ 2j2(-u i '2V il r 1 (i,i') + u* 2 v; 1 r 2 (i,i'))} (9) 
a' 

where 

Bi(i,i') = (—-UiiiVi2sin2a(l — h) + Ui'2Vii(sin 2 a + hcos 2 a)) 

Ai(i } i') = {-^-U* x Vy 2 sin2a{l - h) + [/^^(rfa + /icos 2 a)) 

1 2 1 

r^i') = m Xi m x .,U i2 (-—f 2 {x Xi ,q,x x . l g)V i/1 + £ — A(z', k)T kl f 2 (x Xtik )) 

m k=i m i k 

2 

r 2 (i,0 = ^(-^(^^-^x^^K'i + I] A(i', k)T kl f 2 (x xdk )) 

k=l 

A(i,k) = VnT k i - V i2 T k2 , 
fi{xij) = 1 — — inxy + Inxwj 

X{j 1 

/ 2 (a;,2/) = (— — -/ni V —lfiy) 

x — y x — 1 y — 1 

ar y = m^/mj ; = (10) 
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the various m's represent the masses of the corresponding particle 
U and V represent the matrices which diagonalises the 



m 



ho/™?H - 



In (|), 
and h 

chargino masses. T is the matrix used for diagonalisation of stop mass matrix. 
m is the average mass of the first two generation of squarks. Starting with 
Cq^tuw), there values at a scale m& can be evaluated by solving the RGE 
equations. The result is 



C Qi (m b )=r)-^°C Qi (M w ) 



(11) 



where r\ 



a s (m b ) 



u 



>7q 



-4 and /3 = 11 - (2/3)71/ = 9 | 



a s {M w ) 

Taking the matrix elements of ([?[) between vacuum and \B®) as before, 
we see in effect that Cq 1 and Cq 2 enter multiplied by m Bs unlike C w in (H) 
which entered as miC w . For large tan/3 thus, the decay B s — > u + n~ will 
be dominated by (|7|) and we get : 



-pSUSY 
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The co-efficient C\o in (|T^) has been evaluated in NLO approximation by 
Misiak || in SM.We use his value since SUSY contributions are expected to 
change this value only slightly and also because in the region of large tanfi, 
the RHS of flT2D will be dominated by the Cq 1 ,Cq 2 term. As we can see 
from (p[)- (p] ) the co-efficients Cq 1 and Cq 2 depend, in addition to SM pa- 
rameters^ the MSSM parameters. The most economical version of MSSM is 
the one where the spontaneous breaking of SU(2) U(l) is achieved through 
radiative effects. Further at Unification scale of Mgut{^ 10 16 GeV) scalars 
have a common soft breaking mass term m and Gauginos also have a unified 
soft breaking mass M.Thus at scale Mc^,MSSM in addition of SM parame- 
ters have five parameters : m,M, the ratio of the vacuum expectation values 
of the two Higgs fields tan/3, the trilinear soft breaking scalar term coefficient 
A and the bilinear soft breaking term coefficient B.The values of these pa- 
rameters at scale ~ m b relevant to us can be related to their values at Mqut 
through RGE which has been extensively discussed in literature H,[I2|.As 
pointed out in |TJ| to suppress the large contribution to K° — K° mixing it 
is sufficient to require degeneracy of the soft SUSY breaking mass in squark 
sector. Thus, the strict universality of all scalar masses is not necessarily re- 
quired in context of SUGRA model. We thus consider in our estimate of 
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study B s — > ,in regions of parameter space where the condition of 

universality of scalar masses is relaxed. 

The parameter space of MSSM is constrained by several pieces of known 
experimental results, the W ± , Z one loop mass corrections, (g - 2) of the \i 
meson and most importantly by the lower and upper bounds to b — ► 57 decay 
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B(b -> s 7 ) = (0.6 - 5.4) x 10 



-1 



MSSM parameter space in this context has been extensively analyzed by 
Lopez et.al |g,Goto et.al 0. We present in Fig(|) the rate of T SUSY (B S -> 



/i + /i~) (fL^) relative to their SM values for some of the allowed region of 
the parameter space for large values of tanfi in SUGRA model where the 
universality of scalar masses in SUGRA model is relaxed. The SM branching 
ratio of B s — ► is 2.6 x 10 -9 , so for large values of tanfl the branching 

ratio can be as large as ~ 10 _6 .This is well within the experimental reach of 
the B-factories expected soon. The B s — > fi + fi~ observation or otherwise will 
thus serve as a very useful constraint on the parameter space of MSSM. 

We present the result of SUGRA model in which we are assuming the 
universality of scalr masses in Fig(|]) We also present the result of a more 
predictive class of Supergravity (SUGRA) theories : the Moduli (or No- 
scale) scenario ( m = A = at unification scale) and the dilaton scenario 
(m = M/^/3,A = —1) in Fig(^).We see that with the allowed choice of 
parameters the B s — > n + n~ rate depend crucially on the value of the Higgs 
scalars and pseudo-scalars.For low values of those masses with not too small 
splitting of the stops, the rate can be enhanced by more than two orders 
of magnitude (Fig |3]). Thus B s — > n + n~ data together with estimates of 
masses of superpartners will be very useful in cross-checking these models 
with experimental results. 

We next turn to the decay B s — > l + l~ 7, which was studied in the frame- 



work of QCD corrected effective SM Hamiltonian by Eilam et.al. JET]. As 
is pointed out there , the radiative process unlike the non-radiative chan- 
nel does not face the helicity suppression factor proportional to mi. The net 
result that one ends up is paradoxical in that the Branching ratio for the 
radiative mode turns out to be higher than the one for the non-radiative 
mode notwithstanding the fact that the former has an extra factor of a. 
This is very easily seen by considering the set of graphs where the photon 
line is hooked on to the various external lines of the non-radiative decay pro- 
cess. Graphs when the photon line is hooked on to the lepton gets a factor 
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of mi and so can be ignored. The two graphs with the photon hooked on 
to the b and s quark give rise to the following amplitude for the process 

BsipB.) 



a(b s -> i + r 7) 



+ p R r 

2C 7 m 



Ae, 



tfy- p' s + m s 



-2p s .p 7 
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Q b ] b [C e 9 ff l^l + Cwl 7 M 75 
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P 



[Pr PT 7T 7 Qb 
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+ m s 



(13) 



In flTJD Q s , Q b are the charges of b and s quarks , p b and p s respectively 
represent the momentum of the b and s quark inside the B° s and in the 
constituent quark model can be written as p bs = ps s .The fact that 

' m B s 

m s <C m b allows one to retain only the leading term in (m s /m b ) in eqn(|T3|) 
and gives to the amplitude: 



A(B S -> 1+ r 7) 



e aG F f Bs \ 
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2C 7 m 



7'"B a 



p z 



)(ll,l)+C 10 (l^ l5 l)} (14) 



Where once again we have used eqns(f2|),(|]),((4p. Using the values of C Q ,C 7 
and Cio the Branching ratio into the B° s — > /i + /i~7 mode is estimated at 
4.6 x 10~ 9 [|11|] , almost twice the branching ratio of B® — > ur . 

Let us now examine the effect of the NHB exchange diagram in SUSY 
models on the amplitude A(B S — > / + / _ 7).The contribution can be written in 
terms of Cq 1 and Cq 2 terms. So A(B S — > l + l~ 7) is : 
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NHB , 



B„ 



/+r 7 ) 



otGi 
12V2 



7T 



>>ib 



Vb-V~i 



[c Ql (ii) + c Q2 (h 5 i)) 



(15) 



where we have made the same approximations as in deriving eqn(|TJ]) i.e. 
retaining only the leading power of (m s /m b ). 
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However 



e„(0|S7" ^P R b\B° s ) = \e v f (0\sP R b\B° s ) 

+~ e ^7 (0|^,7l^|5 s °) 
= (16) 

from eqns.([J) and the transversality condition for real photon e u p" = 0. 

Thus,the NHB terms are negligible even for large tan/3 co-efficients Cg and 
Cio do not change much in MSSM. Changes in Cj in SUSY are constrained 
by b — > 57 data.This change too is less than a factor of two. Thus the SM 



estimate of (TA) effectively still holds in MSSM also. 

We thus can have a change in the pattern of conclusion of when the 
Neutral Higgs Boson contribution are included with the large tan/3depending 
upon SUSY and most importantly upon the values of Higgs masses and 
tanj3. The radiative mode B° s — > / + /~7 branching ratio stay put at about 
4.6 x 10 _9 , whereas the non-radiative (suppressed earlier without NHB) can 
jump by a factor of the order ~ 10 2 to a value ~ 10~ 7 .With the flux of 
expected to go up to observe such branching ratios, one thus has a concrete 
and interesting situation to test SUSY-models with large tan(3 in studying 
the radiative and non-radiative dileptonic decays of B® 
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FIGURE CAPTIONS 



Figure 1 : Wilson coefficient Cq 1 as a function of light chargino mass in 
SUGRA model without nonuniversality of scalar masses (pseudo-scalar higgs 
mass is taken to be in range 80 < itia < 700), m = M, A = - 1/2. 
Figure 2 : Wilson coefficient Cq 2 as a function of light chargino masses in 
SUGRA model without nonuniversality of scalar masses with same parame- 
ters as above. Parameters are 80 < uia < 700, m = M, A = — 1/2 
Figure 3 : Ratio of T SUSY /T SM (SUSY rate normalized to SM rate) as a 
function of light chargino mass in SUGRA model without nonuniversality of 
scalar masses. Parameters are 80 < itla < 700, m = M, A = — 1/2 
Figure 4 : Y SUSY /T SM as a function of light chargino mass in 

(a) Dilaton (m = M/VS , A = - 1) 

(b) SUGRA ( m = M , A = - 1/2 ) 

(c) Moduli ( m = A = 0) 

Points with dots (.) are consistent with allowed parameter space whereas 
points marked with plusses (+) are those excluded by our choice of parame- 
ter space (including b — > constraints). 



10 




11 




12 




0.1 







50 



100 



150 



200 



250 



Figure 3: 



13 




14 



